Magnon Supersolid and Anomalous Hysteresis in Spin Dimers on a Triangular Lattice 
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We study the magnetic phase diagram and hysteresis behavior of weakly coupled spin dimers on 
a triangular lattice using the cluster mean-field method with cluster-size scaling. We find that the 
magnetization curve has plateaus at 1/3 and 2/3 of the total magnetization, in which local singlet 
and triplet states form a superlattice pattern. Moreover, if increasing (decreasing) the magnetic 
field from the 1/3 (2/3) plateau, the Bose-Einstein condensation (BEC) of triplons occurs on the 
superlattice background, leading to the transition into magnon supersolid phase. We also find that 
the first-order transition between these solid states and the standard magnon BEC state exhibits 
an anomalous hysteresis upon cycling the magnetic field; the transition can occur only from solid 
to BEC, and the system cannot return to the initial solid state in the reverse process. 

PACS numbers: 75.10.Jm, 75.45. +j, 67.80.kb 



The physics of bosonic elementary excitations in cou- 
pled spin-dimer systems has been attracting considerable 
attention since the observation of magnon Bose-Einstein 
condensation (BEC) in TICUCI3 0t3|- The strong in- 
tradimer antiferromagnetic interaction in spin-dimer sys- 
tems [4H8( provides a quantum disordered (QD) singlet 
state, and the system undergoes a phase transition to a 
magnetically ordered state when applying a sufficiently 
strong magnetic field. This quantum phase transition is 
now well established as a BEC of spin-triplet states with 
S z = +1, called "triplons," moving in the non-magnetic 
singlet background 0, Moreover, some materials ex- 
hibit magnetization plateaus at fractional values of the 
total magnetization where the triplons are crystal- 

lized in a superstructure pattern and form an incompress- 
ible "solid" state. As a next step, it is naturally expected 
that one could observe a more exotic state, namely the 
magnon supersolid (SS) [H-17|, in which the diagonal 
(solid) and off-diagonal (magnon BEC) long-range orders 
coexist fl8j ] . 

The possibility of SS states has been discussed also 
in the contexts of solid helium-4 19l- 22j and ultracold 
Bose gases in optical lattices |23l43d |~ In the former 
case, although Kim and Chan's torsional oscillator ex- 
periments observed superficial nonclassical rotational in- 
ertia in 2004 [Til , it is now interpreted as shear modulus 
stiffening rather than superfluidit y in the latest experi- 
ments with an improved oscillator [22j . In the latter case, 
previous theoretical studies have established that suffi- 
ciently strong dipolc-dipole interactions between atoms 
or molecules, which are long -ranged, allow for the emer- 
gence of SS phases 23l430l|. However, the strength of 
dipolar interactions available with current technology is 
still significantly smaller than temperatures, and SS has 
never been realized yet. Given the presence of a rich 
variety of materials and the large spin-exchange interac- 
tions compared to achievable temperatures, spin-dimer 
systems are advantageous over the other systems for re- 



alizing SS states. 

Nevertheless, no clear SS behavior has been observed 
in any spin-dimer materials so far. This may be be- 
cause the parameter window for SS is very small. Previ- 
ous theoretical analyses with square-lattice models [i"lT - 
[fij l have indeed shown that the magnon SS state can 
emerge only in a very narrow range of parameters (e.g., 




FIG. 1: (color online) Ground-state phase diagram of spin- 
1/2 dimers on the triangular lattice obtained by the CMF-9 
calculations. We set Jo = 2(J2 + J3) = 1. The dashed green 
lines are the phase boundaries determined by the cluster-size 
scaling. Lower panels: schema of the anomalous hysteresis 
process when decreasing and then increasing the magnetic 
field from an initial NS state along the route marked by the 
purple allow in the phase diagram. 
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magnetic field) even when it exists. In this Letter, we 
study the ground-state properties and hysteresis behav- 
iors of weakly coupled spin-dimers on a triangular lat- 
tice. Taking into consideration that in the case of lattice- 
boson systems the SS states on a triangular lattice are 
fairly robust against a domain-wall formation unlike in 
the square- lattice case [24, 2^, 29j| . it is speculated that 
triangular-lattice geometry is more suitable for explor- 
ing SS states also in spin-dimer systems. Furthermore, 
since a novel anomalous hysteresis has been predicted in 
lattice-boson systems (30j . it is imperative to examine 
whether spin-dimers can exhibit such hysteresis. 

We quantitatively determine the ground-state phase 
diagram of the triangular-lattice spin-dimer model even 
far away from the strongly dimcrized limit as shown in 
Fig. [TJ The frustration coming from the competition 
of two interdimer interactions makes numerical analyses 
such as quantum Monte Carlo simulations rather diffi- 
cult and time consuming for frustrated spin-dimer sys- 
tems. Because of this problem, there has been no the- 
oretical (quantitatively reliable) prediction of the com- 
plete phase diagram for frustrated spin-dimer models 
even in the square-lattice case as far as we know [H- 13[. 
To overcome this difficulty, we use a cluster mean-field 
(CMF) method with cluster-size scaling (3l|. As shown 
in Fig. [TJ the magnon SS state emerges in a wide region 
of the phase diagram [in between the two magnetization 
plateaus of the triplon "normal" solid (NS) states]. The 
transition between the solid (NS or SS) and BEC states is 
of first order (thick lines) . We show that the phase tran- 
sition process upon cycling the magnetic field exhibits 
an anomalous magnetic hysteresis (see the schematic il- 
lustration in Fig. [1} analogous to the melting transition 
in a triangular- lattice system of Bosc gases [30[ . In the 
anomalous hysteresis, although the transition from NS to 
BEC can occur (left panel) , in the reverse process the NS 
and SS states are ignored and no reverse transition occurs 
(right panel). This phenomenon is completely different 
from a conventional first-order transition with a hystere- 
sis loop. 

Below we will give a detailed description of the above- 
mentioned physics. The triangular-lattice spin-dimer 
compounds such as Ba3Mn 2 0s p-Q; Sr 3 Cr 2 08 [H, [33j], 
and Ba3Cr208 3J, |35| have a stacked bilayer structure 
shown in the inset of Fig. [T] Assuming here that the cou- 
pling Ji between the dimer planes is weak, we consider 
the following frustrated 5 = 1/2 spin-dimer model (Til — 
13[ on a single triangular-lattice plane: 
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(id) 



-Ja E (*<•' ' %2 + S h 2 ■ S jt i) - H^SU, (1) 



(i,j) 

where Si <a = (S? a ,Sf a ,Sl a ) is the spin-1/2 operator 
attached to the site i, the two spins within a dimer 



are labeled by the subscript a = 1,2, and the second 
and third sums run over nearest-neighbor (NN) sites. 
All the couplings are antiferromagnetic and isotropic in 
spin space; Jo, J 2 , and J3 are the intradimcr, direct 
interdimer, and crossed interdimer interactions, respec- 
tively. The last term is the Zeeman coupling with a 
magnetic field H = gfJ-sh (in energy units). In the 
limit of Jq,H 3> J 2 , J3, the above model is mapped 
onto the hardcore Bose-Hubbard model with NN hop pin g 
t = -(J 2 - J 3 )/2 and NN repulsion V = (J 2 + J 3 )/2 [H, 
which describes hardcore triplons moving in the singlet 
sea. We focus here on the case of t > (J 2 < J3) as 
in usual bosonic models in order to extract the Bosc- 
Hubbard physics in the isotropic spin system given in 
Eq. Q. 

We study the ground-state properties of the spin-dimer 
model in the original form [Eq. (JTJ] for a realistic strength 
of the intradimer interaction, J 2 + J3 = 0.5 Jo, away from 
the strongly dimerized limit. We set Jo = f hereafter. 
Our CMF method [3l[ can take into account the effects 
of quantum and thermal fluctuations by exactly diagonal- 
izing the cluster Hamiltonian with a mean-field boundary 
condition. Here we use the triangular-shaped clusters of 
Ac = I, 3, 6, and 9 dimers (2, 6, 12, and 18 spins) and 
then perform the cluster-size scaling. This method can 
treat quantitatively frustrated systems avoiding the no- 
torious minus-sign problem. 
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FIG. 2: (color online) The magnetization process at J2 — J3 = 
—0.075 with the illustrations of each phase. ^ and rh z are the 
order parameters of BEC and NS phases, respectively. 
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FIG. 3: Examples of the cluster-size scalings of the phase 
boundaries between BEC-SS (triangle, H = 1.7), BEC-NS 
(open circle, H = 1.4; filled circle, H = 2.2), and BEC-QD 
(square, H = 0). 



In Fig. [TJ we already presented the ground-state phase 
diagram in the (| J2 — J3I, i/)-plane obtained by the CMF 
calculation with the nine-site cluster (CMF-9). To ex- 
plain each phase in the phase diagram, we plot in Fig. [2] 
the magnetization curve m z {H) = X^wfi + S*^/M at 
J 2 - J 3 = -0.075 with the BEC (*) and solid \fh z ) or- 
der parameters. Here M denotes the number of dimers. 
At zero magnetic field, the antiferromagnetic intradimcr 
interaction J produces a gapped QD state with m z = 0, 
in which the spins in each dimcr form a singlet state. 
This state cannot be described by a classical spin picture 
(or by the single-spin mean-field approximation). As in- 
creasing the magnetic field H, the energy of the triplet 
excitations with S z = +1 gets lower and the system un- 
dergoes a transition into a magnetic state with closing the 
spin gap at the critical field H c . This phase transition 
can be understood as the BEC of magnetic quasiparti- 
clcs (triplons) @, 0] ■ The order parameter of the magnon 
BEC phase can be written as 'F = 'Yl li {bi)/M using the 
semihardcore boson operator &i = {Sf x — 5 , l + 2 )/v / 2- In the 
original spin language, the BEC of triplons is translated 
as the canted antiferromagnetic state shown schemati- 
cally in Fig. [2j When the magnetic field increases fur- 
ther, the density of triplons gets closer to the commensu- 
rate filling of the triangular lattice, leading to the tran- 
sition into the NS phase. The NS state at the magneti- 
zation plateau with m z = 1/3 has a x v3 solid order 
with the "singlet-singlet-triplet" (s-s-t) superlattice pat- 
tern (Fig. [21 upper middle panel), in which the original 
translational symmetry of the triangular lattice is bro- 
ken. This state is characterized by the order parameter 
m" = + %) exp(zQ • n)/M with Q = (4tt/3, 0). 

Additionally, we can see another plateau at m z = 2/3, in 
which the NS state has the triplet-triplet-singlct (t-t-s) 
pattern (Fig.[5J lower middle panel). At higher fields, the 
system enters into the BEC phase again and eventually 
all the spins are aligned parallel to the field direction. 

In the region sandwiched between the two plateaus, 
both the BEC and solid order parameters can have a fi- 
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FIG. 4: (color online) Anomalous hysteresis behavior at J2 — 
J 3 = —0.101. The solid and dashed lines correspond to the 
(meta)stable and unstable states, respectively. The inset is 
the schematic illustration of the characteristic structure of 
the energy. 



nite value. When applying a magnetic field (or doping 
triplons) on the m z = 1/3 NS state, the BEC of triplons 
occurs in the s-s-t background and the magnon SS state 
emerges. The spins on both sublattice sites cant slightly 
away from the field direction (Fig. [2] upper right panel). 
Similarly, another magnon SS state also emerges by dop- 
ing holes of triplons on the m z = 2/3 NS state. As seen 
in the phase diagram of Fig. [IJ the magnon SS states 
can emerge over a reasonably large region in between the 
two NS states in contrast to the case of the square lat 
ticc [ll 



14[. We also apply a cluster-size scaling of the 



phase boundaries with the scaling factor A (see Ref. [31 
for the definition of A), which has successfully produced 
quantitatively reliable phase diagrams for related hard- 
core boson models [3l|, |3tJ ■ We perform a linear fit of the 
data for the phase boundaries obtained with Nq = 3, 6, 
and 9 (A = 1/3, 1/2, and 5/9), which shows a good fit as 
shown in Fig. [3] The scaled result denoted by the dashed 
green lines in Fig. [1] indicates that the SS region survives 
in the limit of Ac — > 00 (A — > 1). 

Now we present a detailed description of the anoma- 
lous hysteresis behavior already mentioned earlier (Fig.[TJ 
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FIG. 5: (color online) The emergence region of the anomalous 
hysteresis behavior. The left and right panels show the CMF- 
9 and scaled data, respectively. The value of \J% — Jz\ is 
required to lie within the range sandwiched between the two 
vertical dashed lines for the anomalous hysteresis. The green 
dotted curves indicate the metastability limits of BEC phase. 



lower panels). We see in the phase diagram that re- 
entrant (two-step) first-order transitions from BEC to 
solids (SS and NS), and back to BEC can be induced by 
a magnetic field near the tips of the NS phases. However, 
since the actual transition process in first-order phase 
transitions is determined not only by the ground states, 
we need to take into account also the metastable states. 
Figure [3] shows the CMF-9 data of the magnetization 
curve m z (H), including metastable and unstable solu- 
tions, around the m z = 2/3 plateau at J2 — J3 = —0.101. 
The stable (solid curve) and unstable (dashed curve) 
states correspond to the minima and maxima (or sad- 
dle points) of the energy landscape, respectively. When 
starting from an initial NS state, e.g. at point i, and 
decreasing the magnetic field, a transition to BEC oc- 
curs at the metastability limit of SS phase through the 
path i—^a—^x—^x'—^e. In a usual hysteresis, if wc 
perform the reverse process (increase the field strength) 
starting from the BEC state, the system returns to the 
initial NS state through a different path, and a hystere- 
sis loop is formed. However, surprisingly, the BEC state 
keeps the metastability even if increasing H as clearly 
seen in the figure (e — > /). This unidirectional behavior 
is due to the absence of the metastability limit of the 
magnon BEC state. In this hysteresis cycle, obviously, a 
standard hysteresis loop is not formed and the Maxwell 
equal-area rule is not applicable. Although the energy of 
the system cannot be calculated directly in the large-size 
CMF [3l| , the shape of the energy curve can be predicted 
by the magnetization curve m z (H). The curve of the 



energy should form the unconventional structure shown 
schematically in the inset of Fig. [3j in which the line of 
BEC state passes through the loop consisting of NS and 
SS, instead of the usual swallowtail one [H]. Note that 
the same unidirectional behavior can be also obtained by 
increasing first and then decreasing H (through the pass 
i — > b — > y — > y' — > f and then / — > e) and around the 



NS plateau with m z = 1/3. 

The blank (uncolored) regions of SS and NS in Fig. [4] 
are ignored when sweeping the magnetic field from a BEC 
state, leading to the anomalous hysteresis behavior. We 
have to pass through thermal cycling in order to reach 
these solid states. This means, at the same time, that 
the SS and NS states could be missed to observe in ex- 
periments with varying only the field strength at a fixed 
low temperature. For observing the anomalous hystere- 
sis, the initial NS state could be prepared e.g., by cooling 
the sample in a high magnetic field. 

In summary, we have used the cluster mean-field the- 
ory and scaling analysis to study the ground-state phase 
diagram and the hysteresis of the frustrated S = 1/2 
spin-dimcr model on a triangular lattice away from the 
strongly dimerized limit for J3 > J2. The magnetiza- 
tion curve exhibits two plateaus at m z = 1/3 and 2/3, in 
which the singlet and triplet states form a y/3 x \/3 solid 
order with ordering vector Q = (47r/3, 0) (the NS state of 
triplons). We found that more exotic magnon SS states 
emerge in between the two plateaus of the NS states. 
Since the SS phases occupy a broad region of the phase 
diagram in contrast to the square-lattice case, it should 
be easier to experimentally discover these phases in trian- 
gular spin-dimer materials. Moreover, we found that the 
metastability of the system yields an anomalous hystere- 
sis behavior without forming a standard hysteresis-loop 
structure upon cycling the magnetic field. In the anoma- 
lous hysteresis cycle, the actual transition can occur only 
unidirectionally from the solid (SS or NS) phase to the 
uniform magnon BEC phase. The Ginzburg-Landau de- 
scription of the anomalous hysteresis will be reported 
elsewhere 37], 
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